Bottomonium above deconfinement in lattice nonrelativistic QCD 
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We study the temperature dependence of bottomonium for temperatures in the range 0.4T C < 
T < 2.1T C , using nonrelativistic dynamics for the bottom quark and full relativistic lattice QCD 
simulations for Nf = 2 light flavors on a highly anisotropic lattice. We find that the T is insensitive to 
the temperature in this range, while the Xb propagators show a crossover from the exponential decay 
characterizing the hadronic phase to a power-law behaviour consistent with nearly-free dynamics at 
T ~ 2T C . 

PACS numbers: 11.10. Wx, 12.38. Gc, 14.40.Pq 



Introduction - Heavy quark bound states are impor- 
tant probes of the dynamics in the Quark Gluon Plasma: 
charmonium suppression [ll Q has been observed at a 
variety of energies at SPS Qand RHIC Q. While the 
melting of bound states certainly reduces quarkonium 
production, the converse is not necessarily true: differ- 
ent, even competing, effects make it difficult to interpret 
charmonium suppression patterns. It has been noted that 
such effects should be less significant for bottomonium 
(see e.g. Ref. [B[ for a review). Since at LHC energies 
bottomonium will be produced copiously d, 0] , precision 
studies of the suppression pattern and its unambiguous 
link with the spectrum of bound states should be possi- 
ble. The advent of the LHC calls therefore for precision 
studies of bottomonium at high temperature. 

Due to the large mass of the bottom quark, it is cus- 
tomary to study the bottomonium spectrum at zero tem- 
perature using nonrelativistic QCD (NRQCD) fsWloj and 
other effective field theories llj. In this Letter, we em- 
ploy NRQCD to study the response of bottomonium to a 
thermal medium of quarks and gluons in the temperature 
range 0.4T C < T < 2.1T C , at the onset of the initial tem- 
perature attained in heavy ion collisions at the LHC 12 1 . 
We use dynamical anisotropic lattice configurations with 
two light flavors, which have been ex ploi ted before in a 
relativistic study of charmonium 13|, |14| . As explained 
below, the use of NRQCD is a controlled approach which 
avoids many of the unwanted systematic effects encoun- 
tered when using relativistic dynamics for bottomonium 
at nonzero temperature 15l4l7|| . An early study of bot- 
tomonium at nonzero temperature using NRQCD (on a 
quenched background for the 1 So and 3 Si states only) 
can be found in the pioneering work 18j. 

NRQCD at nonzero temperature - In contrast to the 
case at zero temperature [11| , the use of potential mod- 
els to analyse quarkonium at nonzero temperature is less 
well defined due to the uncertainty about which potential 



to use (see e.g. Ref. [19[ and references therein). Recently 
this has been clarified by casting the problem in the lan- 
guage of effective field theory at nonzero temperature 
The series of effective field theories that is ob- 
tained is based on the hierarchy M » T > g 2 M > gT > 
g 4 M, where M is the heavy quark mass and g is the 
gauge coupling. Ref. (2|| provides a clear introduction. 
Integrating out thermal degrees of freedom generates an 
imaginary part for the interquark potential, which high- 
lights the absence of stable states once they are immersed 
in a thermal medium. Limitations of approaches based 
on potential models and the Schrodinger equation are 



discussed in Ref. 27 



In the effective thermal field theory setup |20|42f 
NRQCD is the first theory obtained when integrating out 
ultraviolet degrees of freedom. We study this theory non- 
perturbatively on the lattice and therefore do not require 
weak-coupling arguments as in the hierarchy of effective 
field theories alluded to above. Since NRQCD relies on 
the scale separation M ^> T and we study temperatures 
up to 2T C ~ 400 MeV, its application is fully justified. 

NRQCD has an additional advantage. At nonzero tem- 
perature, spectroscopy for relativistic quarks is hindered 
by the periodicity of the lattice in the temporal direction 
and the reflection symmetry of mesonic correlators, vis- 
ible in, e.g., the standard relation between a correlation 
function and its spectral function, 



G(r) 



doj cosh[w(r - 1/2T)] 
"tt sinh(w/2T) 



p(w) 



(1) 



Nontrivial spectral weight at small w yields a constant r- 
independent contribution to the correlator, which must 
be treated with care [H, [29|. Moreover, it has been 
shown that this contribution can interfere with meson 
spectroscopy [3fj| . which has cast doubt on the status of 
results for the melting or survival of charmonium at high 
temperature 



2 



In NRQCD these problems are not present. Writing 
lj = 2M + u)' and dropping terms that are exponentially 
suppressed when M 3> T [22j], the spectral relation ([T]) 
reduces to 



G(r) 



exp(— bj'T)p{bj') 



(NRQCD), (2) 



2M 



even at nonzero temperature. As a result, all problems 
associated with thermal boundary conditions are absent. 

To study what to expect when quarks are no longer 
bound, consider free quarks in continuum NRQCD with 
energy E p — p 2 /2M. The correlators for the S and P 
waves are then of the form 1221 



Gs(t) 
Gp(t) 



d 3 p 
J2nf 



cxp(-2i? p r) ~ r 



-3/2 



(3) 



(4) 



i.e., they decay as a power for large euclidean time. Of 
course, interactions and finite lattice spacing and volume 
effects are expected to modify this in the realistic case. 

Lattice simulations - Gauge configurations with two 
degenerate dynamical light Wilson-type quark flavors are 
produced on highly anisotropic lattices (£ = a s /a T = 6) 
of size N 3 x N T . A summary of the lattice datasets is 
given in Table I, while more details of the lattice action 
and parameters can be found in Refs. RJlLJ]. We com- 



puted NRQCD propagators on these configurations using 
a mean-field improved action with tree-level coefficients, 
which includes terms up to and including 0(t> 4 ), where v 
is the typical velocity of a bottom quark in bottomonium 
(see Ref. [32| for a discussion of the systematics). The 
states we consider are listed in Table II. 

An accurate determination of bottomonium spec- 
troscopy requires careful tuning of the bare heavy quark 
mass irib to satisfy NRQCD dispersion relations 
Since the main goal of this work is to study the finite- 
temperature modification of NRQCD propagators, an ap- 
proximate choice of a s rrib is made such that nib — 5 GeV. 

Zero temperature results - The zero temperature spec- 
trum from our analysis is summarized in Table [Til This 
spectrum is obtained using a combination of point and 
extended sources in the different channels, in order to 
extract both the ground state and the first excited state. 



N s 


N T 


a' 1 


T(MeV) T/T c No. of Conf. 


12 


80 


7.35GeV 


90 


0.42 


74 


12 


32 


7.06GeV 


221 


1.05 


500 


12 


21 


7.06GeV 


294 


1.40 


500 


12 


16 


7.06GeV 


441 


2.09 


500 



TABLE I: Summary of the lattice data set. The lattice spac- 
ing is set using the IP — IS spin-averaged splitting in char- 
monium 13311. 



state 


a T AE 


Mass (MeV) Exp. (MeV) [34J 




0.118(1) 


9438(7) 


9390.9(2.8) 


2 1 S (Vb(2S)) 0.197(2) 


10009(14) 




l 3 Si(T) 


0.121(1) 


9460* 


9460.30(26) 


2 3 5i(T') 


0.198(2) 


10017(14) 


10023.26(31) 


l x Pi{hb) 


0.178(2) 


9872(14) 




l 3 ^o(xw) 


0.175(4) 


9850(28) 


9859.44(42)(31) 


i 3 A(xw) 


0.176(3) 


9858(21) 


9892.78(26)(31) 




0.182(3) 


9901(21) 


9912.21(26)(31) 



TABLE II: Zero temperature bottomonium spectroscopy from 
NRQCD. The l 3 Si(T) state is used to set the scale. 



Because level splittings are relatively insensitive to m,b 
and to avoid the difficulties in calculating the rest mass 
in NRQCD, we have combined m T (l 3 <Si) = 9460 MeV 
(from the Particle Data Book [34]]) with the mass split- 
tings obtained from our results to predict the bottomo- 
nium spectrum. 

As can be seen from Table [TTJ the zero temperature 
spectrum is reproduced reasonably well. The hyperfine 
splitting between the rjb( 1 So) and T( 3 Si) is much smaller 
than the experimental value, due to the coarse spatial 
lattice spacing, the use of tree-level coefficients, the rel- 
atively heavy sea quarks and contributions of higher or- 
der in v 2 (see Ref. [32j for earlier calculations and dis- 
cussions). In this study we are primarily interested in 
qualitative changes to the T and Xb correlators as the 
temperature changes, so a precision determination is not 
a major concern. 

T and Xb in the plasma - To investigate thermal ef- 
fects, we focus on the T and Xb states, computed with 
point sources. Following the discussion above, our aim is 
to see a transition from exponential decay in the hadronic 
phase, G(t) ~ exp(— AEt), characterizing bound states, 
to power law decay, G(t) ~ r -7 , see Eqs. (01 H]), charac- 
terizing quasi-free behaviour. 

In Fig. 1, standard effective masses, defined by 



m eff (r) = -log[G(r)/G(r-a T )], 



(5) 



are shown for both the T and the Xbi propagators at 
various temperatures. Single exponential decay should 
yield a r-independent plateau. In both cases we find 
that at the lowest temperature, T = 0.42T C , exponential 
behaviour is visible provided one goes to late euclidean 
times. Relevant for the topic of this Letter is that in the 
case of the T, the data at the higher temperatures do not 
show any significant deviation from the low-temperature 
result. On the other hand for the Xbi & strong tempera- 
ture dependence is visible, especially at the two highest 
temperatures, ruling out pure exponential long time de- 
cay. We take this as a first indication that the T is not 
sensitive to the quark-gluon plasma up to T ~ 2T C , while 
the Xb may melt at much lower temperatures. 



3 




1 3 S, (Y) 

x r=o.42r 
o r=i.05r 
□ 7=1. 407 

O 7=2.097. 



40 

xla 



3= 

(D 



) 
: 


i 


i i 

i 3 p, ( Xb1 ) 






x T=0.42T 


8 




7=1.057 


e 




□ 7=1.407 c 






O 7=2.097 c 




e ™ 


I , I 



40 

xla 



FIG. 1: Effective mass plots for the T (above) and Xbi (be- 
low) using point sources for various temperatures: note the 
different temperature dependence. 



To investigate the behavior of the \b propagators in 
more detail, we display in Fig. 2 the Xbo, Xbi and Xbi 
propagators on a log-log scale, at the highest tempera- 
ture. The straight line is a fit of the form G(t) = cr~ d , 
which is motivated by the continuum expression in the 
absence of interactions (U]). We conclude that a power de- 
cay describes the data well at large euclidean time, with 
a power d = 2.605(1) which is close to the continuum 
noninteracting value of 5/2. 

To visualize the approach to quasi-free behaviour in 
another way, we construct effective power plots, using 
the definition 



7eff(T) 



G'(r) 
G(r) 



_G(i 



- G(T - Or) 



2a T G{r) 



(6) 



where the prime denotes the (discretized) derivative. For 
a power decay, G(r) ~ r -7 , this yields a constant result, 
7eff (t) = 7. On the other hand, for an exponential decay, 
G(t) ~ exp(— AEt), this yields a linearly rising result, 
7eff (t) = AEt. The results are shown in Fig. 3. We con- 
firm again that the T displays essentially no temperature 
dependence, while for the Xbi we observe a tendency to 
flatten out, corresponding to power decay at large eu- 
clidean time. Also shown are the effective exponents in 
the continuum noninteracting limit. In the case of the 
Xbi 1 we observe that the effective exponent tends towards 
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FIG. 2: \b propagators on a log-log scale, at the highest 
temperature T = 2.09T C . The straight line is a fit to G(r) = 
CT~ d , with c = 223.2 ± 0.5 and d = 2.605 ± 0.001, using 
r/a T = 10, ...,15. 



the noninteracting result at the highest temperature we 
consider. 

Summary - We have studied the behavior of T and Xb 
at high temperature using anisotropic lattice simulations. 
The bottom quark is treated via NRQCD and the light 
quark dynamics is realized by exploiting lattice configu- 
rations with two flavors of dynamical quarks. The non- 
relativistic approximation for the bottom quark is well 
justified in the range of temperatures we have explored, 
and has many technical advantages. At nonzero tempera- 
ture, the main benefit of using NRQCD over the standard 
rclativistic formulation is that the only temperature de- 
pendence in the NRQCD correlators is due to the ther- 
mal medium, and not due to thermal boundary condi- 
tions. We found that this offers a much cleaner signal for 
the crossover between bound and melted states. It also 
improves the prospects for extracting spectral functions 
inverting Eq. © using the Maximal Entropy Method. 
This is currently in progress. It will also be interest- 
ing to compare the results presented here with those ob- 
tained using a relativistic treatment of bottom quarks, 
employing the same anisotropic action as was used for 
charmonium hi 



This will give an estimate of the 
possible systematic uncertainties inherent in the two for- 
mulations. 

Our results indicate that the T shows no temperature 
dependence up to 2.09T C , while the Xb propagators are 
sensitive to the presence of the thermal medium imme- 
diately above T c . Power- law decay of the Xb propagators 
is visible at T = 1.4T C , while at the highest temperature 
studied, T ~ 2T C , we found consistency with nearly-free 
dynamics. The effective power, defined in Eq. ((6]), is 
temperature dependent and approaches the noninteract- 
ing result at the highest temperature we considered. It 
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FIG. 3: Effective exponents 7eff (t) for the T (above) and xu 
(below), as a function of Euclidean time for various temper- 
atures. The dotted line indicates the noninteracting result in 
the continuum. 



would be interesting to understand the temperature de- 
pendence and crossover between exponential and power 
decay analytically, within the framework of effective held 
theories mentioned above. 
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